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^ 1 The tunnel current of a Luttinger liquid with a finite density of strong impurities is 

I ^^ I ' calculated using an instanton approach. For very low temperatures T or electric fields 

E the (nonlinear) conductivity is of variable range hopping (VRH) type as for weak 

pinning. For higher temperatures or fields the conductivity shows power law behavior 

^ ■ corresponding to a crossover from multi- to single-impurity tunneling. For even higher T 

^^ ' and not too strong pinning there is a second crossover to weak pinning. The determination 

Q>>. ' of the position of the various crossover lines both for strong and weak pinning allows the 

construction of the global regime diagram. 
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o. 

OO , 1. Introduction 

o 

ID electron systems exhibit a number of peculiarities which destroy the familiar 
Fermi-liquid behavior known from higher dimensions. Main reason is the geomet- 
ry ' rical restriction of the motion in ID where electrons cannot avoid each other. As 
j^ ■ a consequence excitations are plasmons similar to s ound waves in solids. The cor- 
responding phase is called a Luttinger liquid (LL) LlMl. Renewed interest in LLs 
arises from progress in manufacturing narrow quantum wires with a few or a single 
conductingchannel. Examples are carbon nanotubes'^, polydiacetylen 1^, quantum 
Hall edges '^ and semiconductor cleave edge quantum wires '^. 

From a theoretical point of view ID quantum wires allow the investigation of 

the interplay of interaction and disorder effects since short range interaction can 

7 
be treated already within a harmonic bosonic theory . Central quantity is the 

interacti on p arameter K which plays the role of a dimensionless conductance of a 

clean LL I^LLJ. The effect of disorder on transport in LLs has been so far considered 

in two limiting cases: „„^,„^^ 

(i) The effect of a single impurity was considered inl^MiULUJ, Here the conduc- 
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tance depends crucially on K. Impurities are irrelevant for attractive {K > 1) and 
strongly relevant for repulsive interaction {K < 1), respectively. For finite voltage V 
and K < 1, the conductance is ~ Vt<~^ ^. These considerations can be extended 
to two impurities. Depending on th e appl ied gate voltage, Coulomb blockade effects 
may give rise to resonant tunneling I^LLiil. 

(ii) In the opposite case of a finite density of weak impurities^ (Gaussian) disor- 
der is a relevant perturbation for K < 3/2 leading to the localization of electrons. For 
weak extern al electric field E the conductivity is highly nonlinear: a{E) ~ ^-c/VE 
| lz | ld | 14 | io | ^^ Jq^ ^^^ finite temperatures T this result goes ove r into the VRH 
expression for the linear conductivity a ~ e"*^ ^^"^ | 1o | 14 | 1o | 1d | 1/ | ^^ higher tem- 
peratures there is a crossover to cr ~ y2-2if lioj 

On the cont rary, much less is known in the case of a finite density of strong 
pinning centers III which we will address in the present paper. In particular we 
determine both the temperature and electric field dependence of the (nonlinear) 
conductivity for this case in a broad temperature and electric field region. The 
main results of the paper are the conductivities dS]) , dS]) , ([7]) and ([8]) as well as the 
crossover behavior summarized in Fig. [TJ 



2. Model and Instantons 

Starting point of our calculation is the action of interacting electrons subject to an 
external uniform electric field E and strong pinning centers. In bosonizcd form the 
action takes the form 



5"=—— / dxdy^idyLp)^ + {dx(p + fx)^ -'^uS{x-x^)cos{2ip + 2kFX,)^ 



'-1 



The phase (/3(x) is related to the electron density p(x) = TT~^{kF+dx'f){l+'2cos{2(p+ 
2kFx)). kp is the Fermi wave vector, r = y/v and / = FK/vH. v and At = hv/T 
denote the plasmon velocity and the thermal de Broglie wave length, respectively. 
The phase field between the impurities can now be easily integrated out leaving 
only its values ip(xi,y) = 4>i{y) at the impurity sites Xi which are assumed to 
be randomly distributed. The action can then be expressed in terms of Fourier 
components (t)i{y) — \t~^ J2uj "/"^Wne"*""^, ^n = 2Tm/\T- Thus 



f^■ \~^ J V^ W„ I \4>i+l,uJn ~ 4'i,uJn\ , l\i |2 , Ii |2\ 4_ ^^n^' 



-f{ai^i+ai)(j)ifi + Ucs dy 1 - cos {2(t)i{y) + 2nat) [ (2) 

where Oj — Xj+i — Xi and ai = kpXi/'K. Since kpa.i 3> 1 below we will assume the 
ai to be random phases but keep the impurity distance Oj approximately constant 
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Next we consider the current resulting from tunneling processes between 
metastable states, assuming strong pinning and weak quantum fluctuations, i.e. 
K <^ 1. The tunneling process starts from a classical metastable configuration 
(pi which minimizes the impurity potential for all values oi y, E = 0. Hence 
(f)i = 7r(ni — ai) where Ui is integer. Among the many metastable states t here 
is one (modulo n) zero field ground state ^° where rii = n^ — X]i<J^'^j-i]g ^^• 
Here Aa^ = ctj+i — aj and [ajc denotes the closest integer to a. A new metastable 
state follows from the ground state by adding integers g^ = ±1 to the n?. 

Next we consider an instanton configuration which connects the original state 
(f)i with the new state (j)i + TT, rii depends in general on y. To be specific, we assume a 
double kink configuration for the instanton at each impurity site: (j)i{y) = (j^i + n ^ for 
\y — yi\ < Di — d, and (f>i{y) = 4>i , for \y — yi\ > Di + d, with a linear interpolation 
between the two values at the kink walls in the regions Ny — y^l — DJ < d. yizt Di 
is the kink/ anti-kink position, d ~ 1/u is the approximate width of the kinks and 
2Di their distance. It is plausible that in the saddle point configuration all yi will 
be the same, an approximation we will use in the following. With Zi = irDi/a the 
instanton action can then be rewritten as 




i+i - Zi) - fa^Zi + s + In 



cosh((z,+i - z,)/2) 2, 

cosh((z,+i + zO/2) 2 

where the sum goes only over impurities with z^ > 0. s is a constant that includes 
the core action of a kink and an anti-kink: s = In(Cau) ^ 1, where \nC/K ^ 1. 

For a given initial metastable state {4>i}i Si is a function of the variational 
parameters {-Di; i = 1,...,N}. The nucleation rate F and hence the current / 

N 

is given by / ex F ex Yi In dDkex:p{—S/h). Here we employ an approximate 

1=0 

treatment in which we assume Di = D = az/ir for k < i < k + m and Di = 
elsewhere, i.e. tunneling is assumed to occur simultaneously through m neighboring 
impurities. The instanton is then a rectangular object with extension ma and 2D 
in X and y direction, respectively. The instanton action can then be written as 



'-'in 



2h 

'k 



z(Jm{k) + ln(l + e ^^) + mis + Intanh z — ) 

V 2 E„J 



(4) 



Here we introduced the dimensionless field strength fa^ /it — E/Ea where Ea — 
l/(Ata^), K = K/nhv denotes the compressibility. (Jm{k) = (j^fe(l) -I- i^fe+m(— 1))/2 
plays the role of a surface tension of the vertical boundaries of the instanton where 
i^k{q) = q^ — 2q(^Aak — [Aafcjc). In the ground state am{k) is equally distributed 
in the interval < a„i{k) < 2 Lt^. The second and the third contribution in ([4]) 
result from the horizontal boundaries of the instanton and include their surface 
tension s/a and their attractive interaction. The last term describes the volume 
contribution resulting from the external field. 

In addition, we have to include a small dissipative term S'bath = ^ 'm-V In z , 
77 ^ 1, in the action in order to allow for energy dissipation L^^. However, we will 



(3) 
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omit ry-dependent terms in all results where they give only small corrections (apart 
from possible pre-exponential factors which we do not consider). 

A necessary condition for tunneling is dSmst/dz < for z — > cx), i.e. (Jmik) < 
mEa/E. The tunneling probability follows from the saddle point value of the instan- 
ton action where z fulfils the condition am{k) — m^ + tanh 2 — 1 + ^ + jj^^ = . 



3. Results and Conclusions 

We discuss now several special cases: (i) For sufficiently large fields E ^ Ea 
the saddle point is Zs « -^ ^ 1 which gives a tunneling probability F oc 
{E/Ea)~^~^e~'^"^'^^^ . The exponent —1 results from the integration around the sad- 
dle point. Because of small K and correspondingly large kink core action, tunneling 
through single impurities {m = 1) is preferred and hence the nonlinear conductivity 
is given by 

a{E) ^ {E/Ea)*'^ e-'^/^, Ea < E < E^,,, = Eae\ (5) 

in agreement with previous results for tunneling through a single weak link '^ if we 
identify e^"'^ ^ t with the hopping amplitude t through the link. The upper field 
strength for the validity of this result can be estimated from DgU = Zgau < 1 since 
the instantons loose then their meaning. Using u — *■ u^s ~ kp{u/kpy^^^^-^'> we 
find Ei_cT ^ {kp/eK.a){u/kFY^'''^~^\ which can be also read off directly from (HJ 
as i?i,cr ~ EaC^ . Classically {K = 0), -Ei,cr corresponds to the case when the field 
energy Ea the electron gains by moving to the next impurity is smaller than the 
pinning energy u/ k. 

At finite temperatures there is a crossover to a temperature T w EaK dependent 
conductivity 

2 g 

a{T,E)^( — ] e-'^/^sinh(— )— , EKa < T < T^,,, ^ Tac' (6) 

when the instanton extension 2Ds reaches At, i.e. for E < EaT/Ta- For tempera- 
tures higher than Ti ^r isolated impurities are weak. Following the arguments of^!^ 
one expects in this region a ^ T^^^^ . 

(ii) In the opposite case of weak fields, E <^ Ea, tunnehng happens simultane- 
ously through many impurities and the saddle point is Zg ^ 1. In this case we can 
estimate the typical surface tension as ak{m) « 1/m for a chosen pair of sites k and 
k + m, respectively ^^. For very large values of m we can treat ni as continuous 
and the saddle point condition gives rris ~ y^Ea/E ^ 1 and Zs ~ sEa/{2E). The 
tunneling probability and hence the current is proportional to 



/ ~ a{E) ^ e-*V^^^, E < Ea. (7) 

If we write the result in the VRH form LiSI / ^^ g-2ma/5ioc ^q ^.^j^ identify the 
localization length ^loc ~ aK/s of the tunneling charges. There is a crossover to 
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a temperature dependent conductivity if At < ZDs, i.e. for E < sTEa/Ta < Ea 
where 

,Ea^ T 



a{T) ^ e-*V^^^^) sinh(^)^ , EKa/s <T< Tjs 

1 Ea 



(8) 



Results (O and ([8|) are in agreement with those obtained for weak pinning! "^1 1 (jji) 
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Fig. 1. Left: Field and temperature dependence of the conductivity in the various regions of the 
T—E plane. Ta, Ti_cr, s, Ea and Ei,cr are explained in the text. The region Ta > T > Ta/s, E < Ea 
is characterized by activated behavior a ~ e"^"' ^sinh(^^)-^ Right: u — T phase diagram of 
the linear conductivity of disordered LLs. For strong pinning, u > u^ ^ kp(kpa)^~^ and T < 
Ti,cr ~ Ta{u/uc)^"'^~^' , Ta / s Separates the VRH from the single impurity hopping regime. For 
T > Ti,cr impurities become weak. For weak pinning, u < Uc, Ti^^ ~ Ta{u/uc)'^''^^~-^^' separates 
VRH from renormalized power law behavior. For T > Ta the power law is unrenormalized. 



If m is not too large (e.g. for large a) we have to take into account the discreteness 
of TO. An instanton solution exists only for m > y^Ea/E. Since Sinst{z{m),m) 
has always a negative derivative with respect to to at ttt. — > ^jEajE + 0, but for 
reasonably large values of s the interval of to with negative derivative is much shorter 
than 1 and hence the optimal hopping length to^ (E) is the smallest integer exceeding 
yjEajE^ which we denote as \y/ EalE^„ . To be more realistic we have to take 
into account the randomness of the impurity distances a^ such that decreasing the 
field (or the temperature), the current jumps by a factor ~ 5-20^/600 _ Clearly, for 
long wires these jumps will average out. 

Finally, we briefly compare the present case of Poissonian strong disorder, 
UeffO ^ 1 with the Gaussian weak disorder, Uefja <C 1 considered in HO io 14 is ] 
In the latter case u and a are sent simultaneously to zero but the quantity 
u^/a ^ ^q"'^ ^ k\ is assumed to be finite, fo denotes the bare correlation length. 
Fluctuations on scales smaller than ^0 rcnormalize ■Co — ^ C ~ k~p^{'^QkF)^^^^~'^^^ ■ ^t 
low T the conductivity is of variable range hopping type ([5]) up to a temperature 
Tloc = ?iw/$ = ra(u/uc)2/(3-2K) ^jjgj.g y^ ^ kpiakp)'^^^- For higher T there is a 
direc t cross over to cr ^ rp2-2K{T) ^j^^gj-g j^ jg j^q.^^ renormalized by disorder fluctua- 
tions I ° l ^ i This renormalization disappears only at much higher Ta ~ Hv/a. Both 
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weak and strong pinning theories should roughly coincide for u — > Uc ~ kp{akp)^~^ 
where Ta ~ Tier ~ Tioc which is indeed the case since ^ « a. In the strong pinning 
region ^ continues as ^ ~ a/s. 

Experimentally, a linear variable range hopping conductivity has been seen in 
carbon- nanotubes 1^ and polydiacetylen — . 
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